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Abstract. The observed molecular properties of a sample of FIR- luminous and OH megamaser 
(OH-MM) galaxies have been investigated. The ratio of high and low-density tracer lines is found 
to be determined by the progression of the star formation in the system. The HCO + /HCN and 
|. HCO + /HNC line ratios are good proxies for the density of the gas, and PDR and XDR sources 
D ■ can be distinguished using the HNC/HCN line ratio. The properties of the OH-MM sources in 
the sample can be explained by PDR chemistry in gas with densities higher than 10 5 ' 5 cm -3 , 
I , confirming the classical OH-MM model of IR pumped amplification with (variable) low gains. 
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^ | 1. Introduction 

High-density molecular gas plays an important role in the physics of (Ultra-)Luminous 
Infrared Galaxies ((U)LIRGs), giving rise to spectacular starbursts and possibly provid- 
ing the fuel for an active galactic nucleus ( AGN) . The emission lines emanating from the 
7-H ■ nuclear gas provide information about the physical properties of the nuclear environment 
in these systems, e.g. the (column) density, temperature and chemical composition of the 
gas, and the type and strength of the central energy source. It can also provide us an 
insight into the processes influencing the (gas in the) nuclei: the star formation rate and 



£f~) • hi story, fuellin g of a p ossible central black hole and feedback processes. 



Baan et al. I (120071 ) present data of the CO, HCN, HNC, HCO+, CN, and CS line 
emissions of a representative group of 37 FIR-luminous and OH megamaser (OH-MM) 
' galaxies and 80 additional sources taken from the literature. In this work, the molecular 
characteristics of this sample are explained using several models. First, the properties of 
the different ( density) componen ts of the nuclear gas are explained in terms of starburst 



evolution (see iBaan et al. 1 120071) . Als o the chemical properties o f the high-density gas 
' are ana lyzed, using chemical ne tworks (jMeijerink fc Spaans 2005) and radiative transfer 



models iMeijerink et jjD (|2007h . Here only the J=l-0 transitions of the molecules ar e 



considered. A more detailed analysis will be presented in Loenen. Baan fc Spaans I (|2007h 



2. Starburst evolution 

The gas in galaxies is build up from multiple components, each with different densities 
and temperatures. From our data we can derive information about these different com- 
ponents. The CO(l-O) line traces the large scale low-density (critical density 10 s 
cm~ 3 ) component (LD), whereas the lines of HCN, HNC, and HCO + (all n cr >10 5 cm~ 3 ) 
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Figure 1. left, from top 
to bottom: Integrated line 
ratios HCN(1-0)/CO(1-0), 
HNC(1-0)/CO(1-0), and 
HCO+(1-0)/CO(1-0) versus 
FIR luminosity. Squares 
represent reliable values and 
triangles upper or lower 
limits. Filled symbols are 
sources with OH-MM activ- 
ity. 

right, from top to bottom: 
Three FIR luminosity curves 
used in the simulations with 
different maximum lumi- 
nosities, the luminosity of 
a high-density component, 
and the variation of the 
high- versus low-density 
ratio during the outburst 
starting at the upper part of 
the curve. 
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trace the high-density component (HD) which mostly resides in the cores of the galaxies, 
at the sites of the star-formation activity. 

On the left hand side of Fig. [T] the relative contributions of these two density compo- 
nents are shown, by plotting the ratios of the high-density and CO(l-O) line strengths. 
The figure shows that the distribution of line ratios for all molecules increases with FIR 
luminosity, which gives an upwardly curved lower boundary for the distribution at higher 
Lpm- The highest values are found at Lpm^lO 10 ' 5 L Q . The figure also shows that in 
general the OH-MM sources have a much larger spread in HD/CO(1-0). 

This behavior can be explained as the result of ongoing star formation. The FIR 
luminosity of the ULIRG during the evolution of the outburst reflects energy generated 
by the star formation activity. The FIR luminosity integrated over the course of the 
outburst would reflect the amount of high-density molecular material consumed by the 
star formation process and destroyed or r emoved by feedbac k. 

In the following simplified scenario (see lBaan et al. IboOTt ). we employ a model for the 
time-evolution of the high-density components in a galaxy during a FIR outburst. In the 
absence of a representative FIR light curve, we use a diffusion-like expression in time t 
as a response to a starburst starting at t=0 defined as: 

, x 2.5 

L Fm (t) = 1.35L FIR (0) f-j e- ( / T , (2.1) 

where Lfir{Q) is the maximum luminosity of the burst and T is the timescale of the 
outburst. We note that a diffusion curve may not be the most appropriate representation 
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but this curve does resembl e the outcome of starburst-driven FIR evolution 

l2006h . 



of L 

simulations (|Loenen, Baan fc Spaans 



The high-density component HD can be defined as [3 LD, the low-density component. 
As a result, the high-low-density ratio varies with time during the FIR outburst as: 



HD(t) 
LD 



1 
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/ L Fm (t)dt 



IR.int 



(2.2) 



where 7 is the fraction of the initial HD component that is consumed during the whole 
outburst, and LnR.int the FIR luminosity integrated over the whole course of the out- 
burst. The large-scale low-density component LD is assumed to remain unchanged. 

The results of these simulations have been presented on the right hand side in Fig. [TJ 
The top panel shows the FIR light curve of the outburst for three peak luminosities. 
The middle panel shows the luminosity of a representative high-density component for 
the three L-pm, curves. The bottom panel shows the high-low-density ratio for these same 
FIR light curves. Combining the results of this simulation with the data shows that the 
HD/CO(1-0) data points in the left panels are a measure of the evolution of the starburst. 
This implies that the OH MM sources are galaxies in an early stage of star formation, 
which is consistent with OH MM sources b eing found in starburst-dominated galaxies 
( Genzel et aDll998t ISaan fc Klockner 112007t h 



3. Chemistry 

The model presented in the previous section does well in explaining the evolution of 
the different gas components, but it makes no distinction between the different high- 
density tracers. Even though the emission of the different molecules originates in the 
same regions, the line strengths are influenced by the environmental properties like the 
(column) density, temperature, and the type and strength of the prevailing radiation 
field. In order to study the effects of these parameters on the emission characteristics of 
the sources, we remove the intrinsic difference in line strength between the galaxies in 
our sample and use line ratios to find diagnostic properties. Fig. [2] presents the ratios of 
the integrated lines of HCO+/HCN, HNC/HCO+ and HNC/HCN against each other. 

In order to interpret the behavior of the sources in this diagram, we compare the data to 
the theoretical models that treat the chemistry and radiative transfe r of molecular clouds, 
including all the relevant he ating, cooling and chemical processes (see lMeiierink fc Spaans 



2005t iMeiierink et a/.ll2007l, and Spaans e t al. in these proceedings) . A large grid of mod 



els was created bv lMeijerink et all ( 2007 ). varying the strength of the radiation field, its 



type (UV and X-ray), the gas density and the column density. These results are compared 
to our observational data in Fig. [5] (note: not all models are shown, some fall out of the 
range of our figure). 

3.1. PDR models 

The results of the PDR models (UV radiation field) are shown in Fig. [5]with heavy lines, 
where the line styles indicate different gas densities (solid: n—lQ 4 - 5 , dashed: rt=10 5 , and 
striped: n=10 5 ' 5 cm -3 ). The tracks vary as a function of column density, which ranges 
from 7V=10 22 cm -2 (the column density below which the strength of the emission lines 
decreases rapidly) to 7V=10 23 , N=10 23 5 and 7V=10 24 cm -2 for the ?i=10 4 - 5 , ra=10 50 , 
and n=10 5 5 cm~ 3 models, respectively (corresponding to a cloud size of 1 pc; indicated 
by the symbol at one end of the tracks) . Two different radiation field strengths are shown: 
F\jy=1.6 erg cm -2 s _1 , indicated by a plus symbol at the highest column density point, 
and i*uv = 160 erg cm~ 2 s^ 1 , indicated with a circle. 



4 



Loenen, Baan & Spaans 



1 1 1 1 1 1 1 


® J> 


i 

c 


1 1 


□ 

: — 1 — 1 — 1 — 


« « ♦ *, 

■4 \,,, 
1 1 1 h 


- ' * 5- ^ 

r* » r+ 

(pi 

1 1 1 




f ■ 

i i 


9---' 
• • " + + . 


□ 


-0.6 -0.4 -0.2 0.0 0.2 0.4 

logHNC(l-0)/HCN(l-0) 


- ® & D u a n « 
©□ a 

■ 

i i 


► 

i i i i 



-0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 



log HNC(1-0)/HCO + (1-0) 



Figure 2. top left Integrated HCO + (1-0)/HCN(1-0) versus HNC(1-0)/HCO + (1-0) ratios, top 
right Integrated HCO + (1-0)/HCN(1-0) versus HNC(l-0)/HCN(l-0) ratios, bottom left Inte- 
grated HNC(l-0)/HCN(l-0) versus HNC(1-0)/HCO+(1-0) ratios. Explanation about the plot 
symbols and line styles is provided in Sect. [3] 

Two observations can be made, when comparing the data and the models. First of all 
one can see that the models are separated based on the density (n) in the HCO + /HCN 
and HCO+/HNC line ratios. The HNC/HCN line ratio shows no differentiation. This 
can be explained in terms of the critical density of the individual molecules. HCO + has 
a critical density of about 3xl0 5 cm -3 , whereas the critical density of HCN and HNC 
is around 3xl0 6 cm -3 . Therefore the excitation of HCO + will differ from HCN and 
HNC for different densities. A second observation that can be made is that most of the 
OH MM sources cluster together in an area traced by PDR models that have a high- 
density (n^lO cm" 3 ), and a high column density TV^IO 22 cm -2 . This points to the 
classical OH MM model of IR (UV radiation reproc essed by the surrounding gas and 
dust) pumped, low (and variable) gain amplification (Baan 1989). 

3.2. XDR models 

The results of two XDR simulations are also shown in Fig. [31 using thin lines. Again, the 
line styles indicate different gas densities (solid: n=10 5 ' 5 , and dashed: n.=10 6 cm" 3 ); 
column densities range from N—10 22 to 7V=10 24 -10 24 5 ; and the radiation field strengths 
are Fx=1.6 (plus), and Fx=160 erg cm" 2 s" 1 (circle). 

The XDR models are not as well differentiated as the PDR models. Because X-ray 
photons penetrate the molecular cloud much easier, the XDRs do not show the strong 
density dependency seen for the PDRs, making the distinction between different XDR 
models very difficult. The addition of higher transitions and other molecules (e.g. CN, 
CO+, HOC+) will most likely break this de generacy. Another problem with trying to 
identify XDR sources is that they are in general smaller that PDR sources and thus are 
more affected by beam dilution, especially in single dish observations like ours. This will 
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affect the detection rate of XDR sources and could make hybrid sources look like PDR 
sources, even if the XDRs are energetically more important than the PDRs. 

Despite these drawbacks, the XDR models do provide diagnostics, since they clearly 
separate from the PDR models in their HNC/HCN line ratio. The PDR models approach, 
but never cross the HNC/HCN=1 line and the XDR models all have HNC/HCN>1. 
Therefore, the HNC/HCN line ratio is an excellent way to determine whether a system's 
chemistry is dominated by UV or X-ray photons. 

3.3. Terra Incognita 

Not all our observational data in Fig. [2] is covered by the models. The few sources in 
our sample with known H2O MM activity (indicated by plot symbols surrounded by 
circles) are also located in this area, which is characterized by lower HNC and higher 
HCO + line strengths compared to HCN. The fact that this area is not covered by the 
models suggests that other processes influence the line ratios, such as strong shocks. 
Shocks are not treated in the models and can have pr ofound effects on the chemistry in 



molecular clouds as they may selectively destro y HNC ( Schilke et al. II1992I ) and enhance 



HCO+ relative to HCN (jPickinson et al. Ill980h . A simultaneous increase in HCO + and 



a decrease in HNC shifts the PDR models to the uncovered region, implying that the 
H2O MM sources in our sample are UV driven systems with strong shocks. This would 
suggest that these H2O MM sources are similar to shock-induced Galactic H2O maser 
spots (see other contributions in these proceedings). 



4. Conclusions 

Molecular line emissions of multiple species (and transitions) provide excellent diag- 
nostics for understanding the status of the nuclear gas in extra-galactic sources. The 
different molecules trace different (density) components and the ratio of high and low- 
density tracer lines follows the star formation activity in the system. Comparing different 
high-density tracers tells a lot about physical characteristics of the gas. The HCO + /HCN 
and HCO+/HNC line ratios are good proxies for the density of the gas, due to the differ- 
ent critical densities of the species. PDR and XDR sources can be distinguished using the 
HNC/HCN line ratio: PDR sources all have ratios lower than unity and XDRs have ra- 
tios larger than 1 . OH MM sources cluster in a particular part of the diagnostic diagram, 
which is only traced by PDR models with densities higher than 10 5 ' 5 cm" 3 , confirming 
the classical OH MM model of FIR (UV radiation reprocessed by the surrounding gas 
and dust) pumped amplification with low but variable gains. 
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